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Abstract Cell wall directions of tracheids in transverse sec¬ 
tions of Agathis sp. and Cryptomeria japonica and fiber 
bodies in Magnolia obovata were quantified by Fourier 
transform image analysis (FTIA). Three kinds of informa¬ 
tion about cell shapes were determined from the angular 
distribution functions of the net map-power spectrum pat¬ 
tern (NM-PSP): (1) the peak angle expresses the most fre¬ 
quent cell wall direction, (2) the area of the peak is 
indicative of cell wall length, and (3) the dispersion of the 
peak indicates the fluctuation of cell wall orientation. Fluc¬ 
tuations were expressed as the standard deviation (a) by the 
use of statistical analysis. Cell wall directions in Agathis , 
that is, directions of the tangential and radial walls (three 
orientations were detected) were 0° (o = 7°), 59° (12°), 91° 
(9°), and 122° (12°), respectively. Pentagonal and hexagonal 
models were reconstructed as cell shape models in a calcu¬ 
lated ratio of 60:40. On the other hand, directions were 0° 
(5°), 76° (12°), 89° (5°), and 102° (12°) in a calculated ratio 
of 45:55 in Cryptomeria. Fiber bodies were extracted in 
Magnolia and their wall directions were 0° (5°), 50° (7°), 90° 
(6°), and 131° (7°), respectively. The ordinary accumulated 
image analysis (OAIA) was checked at the same time in 
Agathis , and indicated that FTIA could be used as a method 
of reconstruction of cell shape. 

Key words Tracheid shape • Fiber shape • Fourier trans¬ 
form • Cell wall directions • Cell shapes 


Introduction 

The most important structure of wood cells in softwood is 
the tracheid and the fiber in hardwood as ground tissues, 
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although there are various wood cells. Each represents an 
important characteristic of the wood. Until now the shapes 
of these cells, however, were not measured precisely. Fujita 
et al. developed the Fourier transform image analysis 
(FTIA) procedure, consisting of a pre-fast Fourier trans¬ 
form (FFT), which clarifies the shape factor by the dot - 
map and net - map methods, the autocorrelation function 
from the power spectrum pattern (PSP), 2 and the process of 
post-FFT, such as polar coordinate analysis. Then they 3 
detected the cell arrangement and measured the direction 
of cell walls (cell sides) according to the process, and tried 
to reconstruct the most frequent (abundant) cell shape from 
the PSP maximum power. 2 The accuracy and verification of 
analysis are unknown. The controlling process of ordinary 
accumulation image analysis (OAIA) was examined 4 
to prevent various errors and quantitative analysis in FTIA 
was introduced for highly precise quantitative evaluation. 5 
We precisely examined cell arrangements of Agathis sp., 
Cryptomeria japonica , and Magnolia obovata , and deter¬ 
mined the grade of fluctuation of it under this process. 6 Dot 
models and block models were restructured as part of the 
process. In this report, the most frequent angle of the cell 
walls (cell sides) was acquired and reconstruction of the cell 
shapes in transverse sections was performed using the same 
sample images. The grade of fluctuation was also examined 
and both data from FTIA and OAIA in Agathis were com¬ 
pared in detail. 

A number of terms are not defined in this report because 
of their definition in previous reports. 4-6 Terms specifically 
defined in this article are shown in Table 1. 


Materials and methods 

Specimens and picture sampling 

Agathis , Cryptomeria , and Magnolia were analyzed accord¬ 
ing to our previous report. 6 The same sections and input 
methods 6 were used in the present study except that micro¬ 
scopic magnification was higher (10X). Pictures of Agathis 
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Table 1. Abbreviations and definitions of terms 


Abbreviation Definition 


^ol 

ER oi 

P/c=ibw(0) 
Tt=mBw(0) 
P/c=hB\v(0) 
Pac(<P) 
r i> r f , r 2 

A 

- r *-peak 

Lr x , Lr f , Lr 2 


Overlapping error 
Error ratio by overlapping 

Angular distribution function along the low-frequency band width (17 ^ k ^ 64) 

Angular distribution function along the middle-frequency band width (65 ^ k tk 256) 

Angular distribution function along the high-frequency band width (257 ^ k kk 512) 

Angular distribution function after angle compensation 

Type of radial walls on cell shape models. Radial wall 1 and floor-type radial wall and radial wall 2, respectively 
Peak area. Ar 1? Ar f and Ar 2 are peak areas of r l9 r f , r 2 , respectively 
Lengths of r 1? r f , r 2 walls 


and Cryptomeria were directly input to Luzex III (Nireco) 
through a monochrome television camera (Ikegami 
ITC530), to give a one-pixel width on the CRT screen of 
1.39//m, with an input area of 1.42 X 1.42 mm. There were 
about 35 tracheids along the radial direction, and mainly 
contained an earlywood zone and a small area of the late- 
wood zone in Agathis. Latewood-like zones were observed 
repeatedly, although Agathis is one of the tropical trees. 
One thousand tracheids were sampled at a time. There were 
about 20-25 tracheids along the radial direction and the 
total number of cells was about 700 in Cryptomeria. They 
consisted of almost all earlywood cells. Magnolia was 
input to Luzex III through a film scanner Polascan 4000 
(Poraroid, 1280 dpi) and image processing software, 
Photoshop.v.5.0J with a pixel width of 0.83 /um (sample area 
0.85 X 0.85 mm). On the other hand, for the OAIA examina¬ 
tion of Agathis sections, only the earlywood area was se¬ 
lected and input under the same conditions as Magnolia. 


Extractions of the structure element and shape factor, and 
their PSP analysis 

A transverse image of Agathis was transferred to a binary 
image as a wall map (WM) 4 (Fig. la), and simplified to a net 
map (NM) 4 (Fig. lb) after removal of rays according to 
previous reports. 2 4 The same treatments were performed 
for Cryptomeria. Furthermore, the NM was divided into 
fragment maps (FM) 4 of tangential walls (Fig. 2a) and of 
radial walls (Fig. 2c) by cutting off at the corner map (CM) 4 
in order to measure the length and direction of cell sides by 
OAIA. A selective net map (sNM) (Fig. le) was created for 
Magnolia (Fig. Id) as previously described 6 for extraction of 
only the body part of the wood fibers. 

These maps were transformed to the PSPs by FFT pro¬ 
cessing (Fig. lc,f). These PSPs were used for the polar coor¬ 
dinate analysis as programmed by Maekawa et al. 2,3 and 
angular distribution functions were obtained (Figs. 3, 4, 5a). 


Results and discussion 






Fig. la-f. General procedure of Fourier transform image analysis 
(FTIA) for direction analysis, a Wall map (WM), b net map (NM), 
c power spectrum pattern (PSP) of b in case of Agathis sp. d WM, 
e selective NM (sNM), called a body-NM, and f PSP of e in case of 
Magnolia obovata. Note that NM-PSP (c) provides spreading power 
from the origin to the high-frequency bandwidth in a star-like fashion 
in eight directions in Agathis. On the other hand, sNM-PSP (f) shows 
that the power is distributed in a circular fashion and some spots are 
discerned on the circle (arrowheads) and further spots exist on the 
higher frequency power ring (arrows) in Magnolia 


Overall information of NM-PSP 

Structure elements on a NM are cell sides and they contain 
three types of shape factors, that is, length, direction, and 
arrangement. When NM-PSP of Agathis and Cryptomeria 


are surveyed in the low-frequency band width (IBW, 17 ^ k 
^ 64), there are signs similar to DM-PSP on which arran¬ 
gement information appears. We can see star-like rays 
spreading in eight directions from the origin to the middle- 
frequency and high-frequency band widths (mBW: 65 ^ k 
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Fig. 2a-d. Fragment maps and PSPs of Agathis sp. a Fragment map of 
tangential walls (t-FM), b PSP from a, c fragment map of radial walls 
(r-FM), and d PSP from c. Note that power concentrates on the equa¬ 
torial direction in b. On the other hand, power fans out from the 
meridian direction in d. The direction of power crosses the direction of 
cell sides perpendicularly 


^ 256, hBW:257 ^ k ^ 512). These power radiations in¬ 
form us of the length (/) and direction ( 0 ) of the wall sides 
(Fig. lc). On the other hand, eight spots are observed along 
a ring in the mBW on the PSP of Magnolia sNM (Fig. If). 
The ring corresponds to the paired lines of five-pixel width 
and spots on the ring originate from the wall direction (0). 

Concentration and diffusion of the power were calcu¬ 
lated more quantitatively, that is, P k=hBW (0) in Agathis and 
Cryptomeria , and 7\ =mBW _ hBW (0) i n Magnolia (Figs. 3, 4, 5a), 
in order to examine only cell side lengths and directions 
without arrangements. 

Here, we note that three kinds of keys relating to cell 
shape are expressed on the peaks of F(0). Namely, key 1 is 
0 max . It shows the most frequent direction (0 max ) of the cell 
wall, i.e., 0 max = 0 max ± 90°. Key 2 is the peak area. It has a 
proportional relation to the integrated length of each cell 
wall, although the proportionality coefficient is not known. 
Because a cell side that has a one-pixel width is treated as a 


Angular distribution function of NM-PSP 
Agathis sp. P ae 25 , SiS5 , 2 ( 0) 
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Fig. 3a-c. The NM-PSP analysis of Agathis sp. a Raw angular distribu¬ 
tion function, b angular distribution function after angle compensation, 
and c angular distribution after peak separation composed of five 
Gaussian peaks. Compensation was done by using the test patterns 
similar to the structure element of NM (b). It was used in fluctuation 
analysis because c consisted of normal distribution curves. R (correla¬ 
tion coefficient) = 0.997 

o function and there is no power decrease, that is, no miss¬ 
ing order, the relation is kept in whole bandwidth. Key 3 is 
the peak dispersion expressing the fluctuation of wall direc¬ 
tions from 0 max . 
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Angular distribution function of NM-PSP 
Magnolia obovata P 512 ( 0) 


Angular distribution function of NM-PSP 
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Angular distribution function of NM-PSP 
Magnolia obovata P ac , 29S * a5 , 2 ( 
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Angle[° ] 

Fig. 4a-c. The NM-PSP analysis of Magnolia obovata. a Raw angular 
distribution function, b angular distribution function after angle com¬ 
pensation, and c peak separation composed of five Gaussian peaks. 
Compensation was done by using the test patterns similar to the struc¬ 
ture element of NM (b). Even power indicates that there are many cell 
sides not identified in c . R = 0.978 
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Fig. 5a-c. The NM-PSP analysis and reconstruction of Cryptomeria 
japonica tracheid cross shape, a Angular distribution after peak separa¬ 
tion composed of five Gaussian peaks, b pentagonal (house) model, 
and c hexagonal model. Cell side directions and standard deviations o 
are added. Values in parentheses are calculated data from models. 
House model: hexagonal model = 45.2: 54.8 


Experimental conditions for transverse shape analysis 


The FTIA enables us to pull out certain features from am¬ 
biguous structures when many samples are prepared. We 
tried to input as many tracheids as possible in order to 
enhance the reliability of the angular distribution function. 
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For this purpose, lower magnification of picture sampling is 
desirable within the optical resolution R 0 and discrete reso¬ 
lution R d . 4 A wider sampling area sometimes contains un¬ 
desirable heterogeneous parts. Therefore, in this report, we 
decided not to make the observations under such low 
magnification and repeated the samplings mainly from the 
homogeneous earlywood zone. As the PSP from each NM is 
very weak, low power is omitted by digitally processing of 
the gray levels. This problem was overcome 5 by overlapping 
different NMs from several areas, although some convolu¬ 
tion techniques were introduced. 6 The overlapped NM 
image will generally become very complicated. This 
procedure, however, is expected not to bring serious 
influences, being based on the Fourier transform theory, in 
which phase information of the spatial domain is eliminated 
in the power spectrum. Fujita et al. 7 verified this character¬ 
istic by the method of optical Fourier transform, and re¬ 
ported that overlapped sections also generated the original 
periodicity in diffraction. This characteristic was also 
applied to electron diffraction of cellulose crystals. 8 
Midorikawa and Fujita 9 cut off maps of only the earlywood 
zone of softwood into long and slender strips along the 
tangential direction and overlapped them. Then they dem¬ 
onstrated that they could increase the number of tracheids, 
holding the uniformity of a sample. The negative influence 
of overlapping (E ol , overlapping error) must be estimated 
beforehand. 

This error must originate from the overlapped points 
between cell sides on two NM sheets (NM l and NM 2 ), be¬ 
cause these points having two times the contrast are again 
transferred to binary data. When the ratio of lines are m 1 
and m 2 on NM X and NM 2 , the probability of overlapping is 

r\ 

m 1 X ra 2 , namely m in the homogeneous NMs. The line 
position is only a few percent in the area on each NM. The 
overlapping probability of over two sheets, that is, m over2 
becomes negligible when this theory is expanded to numer¬ 
ous sheets (NM„) and the line ratio m is under 0.1. The error 
ratio of overlapping ( ER ol{l+2+3+ +n) = ER oln ) is determined 
by the following formula. 

ER otn = (n - l)m/2 (1) 

The curves of ER ol in the case of m = 0.01, 0.05, and 0.1 are 
shown in Fig. 6. In this report, m was about 0.05. The error 
ratio does not exceed 0.1 even if five sampling images are 
overlapped. 

The most severe problem in cell shape analysis is the 
angular error that is caused by oblique lines. 5 Power is un¬ 
derestimated by around 50% at oblique directions, which is 
unavoidable in the present FFT systems. The degree of 
underestimation of powers was detected by using test pat¬ 
terns similar to cell side fragments that are the structure 
element of NM in order to minimize the negative influence. 
This modified angular distribution function that performed 
angle compensation was written as P ac (0). One thousand 
line segments with 1-pixel width and 30-pixel length were 
spread randomly on a screen as the test (Fig. 7a). These 
segments can be treated as a a function. 4,5 From the PSP of 
this map (Fig. 7b), angular distribution functions were ob- 



Error ratio by overlapping 
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Number of net maps 


Fig. 6. Examination of error ratio caused by overlapping ( ER ol ). In the 
case of Agathis sp. that has an area ratio of about 0.05, five maps can be 
overlapped at less than 0.1 error ratio 
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Fig. 7a-c. Examination of angle error by using fragment test patterns, 
a Many fragments of constant length are dispersed randomly, b PSP of 
a, and c angular distribution functions of b along the middle-frequency 
and high-frequency bandwidths. The y-axis expresses the ratio of 
power level considering the level at 0° as 1 in c. Curves on the high- 
frequency bandwidth were used to compensate for the raw angle distri¬ 
bution functions because the raw functions were traced along the same 
bandwidth 


tained from the same frequency bandwidth as that of the 
real NM-PSP. This operation was repeated 50 times (a total 
of 50000 lines). We compensated the original angular distri¬ 
bution function P(0) to the modified function P ac (0) using 
the curve (Fig. 7c). Because the PSP outputted power with 
point symmetry, the 0^0^ 180° display was sufficient 
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for P ac (0). Peaks were detected at around 0° and 180° in 
both species, and Figs. 4, and 5 show the range of —30° ^ 
0 ^ 210 °. 


Direction analysis of cell wall and reconstruction of 
transverse cell shape 

Reconstruction of tracheid shape in Agathis was performed 
from 0 max values of P ac (0) using key 1. The wave separation 
software Sigma Plot, v.5.0 was used on the assumption that 
each peak has a normal distribution curve in this analysis. 
Maximum directions of peak (0 max ) after the separation 
were slightly different from those before the separation in 
Agathis (Fig. 3). They were adopted for the model recon¬ 
struction. The tracheid arrangement models (block models) 
that were previously proposed 6 were used as the base of 
reconstruction here and we added the information of cell 
wall directions (0) to them by the 0 max . The FM and FM-PSP 
of Agathis are shown in Fig. 2 and demonstrate the relation 
between cell wall direction (6) and power direction ( 0 ). 
When FM is compared with FM-PSP, power is concentrated 
on the direction (0) (Fig. 2b,d) perpendicular to the direc¬ 
tion of cell sides ( 0 ) (Fig. 2a,c). It is clear that NM-PSP (Fig. 
lc) consists of t-FM-PSP (Fig. 2b) and also r-FM-PSP (Fig. 
2d). Here, it must be remembered that the radial direction 
was set by the horizontal direction of the screen (0 = 0°) 
and 0 = 0 ± 90° in the frequency domain (Fig. 2b,d). 

Tangential walls were set around 0 = 90° in the picture 
sampling. It is clear that sharp peaks at 0 = 0°, 180° are 
generated from tangential walls shown in Fig. 3. On the 
contrary, as for radial walls, six radiations were observed 
(Fig. lc) and three peaks at 59°, 91°, and 122°, and their 
point symmetric peaks could also be identified (Fig. 3c). 
Direct interpretation of this NN-PSP is an octagonal cell 
shape. Such cells are rare on the actual transverse section, 
so we can conclude pentagonal cells in the shape of a house 
composed of two walls (tangential walls), one flat floor (r f ) 
corresponding to 0 = 91° (Fig. 3c), and two oblique roofs of 
0 = 59° (r 2 ) and 122° (r x ) (Fig. 8a). The same method was 
used in the case of Cryptomeria. The data are shown in 
Table 2. 

The floor is not of a strict geometry but contains various 
cases. At the tracheid tips, radial walls are often flat and 
their tangential walls are very short. Transverse cell shapes 
are a slender rectangle. We omit such cases in the discussion 
because of their rare occurrence. Tracheids neighboring 
rays have a straight radial wall, and result in the house-like 
shape. When radial files are in contact with each other, 
many tracheids have a typical roof shape, composed of two 
sides. The two sides, however, were often uneven in length 
and the longer ones became flat (Fig. 8b). Moreover, the 
included angles between the two radial walls are often very 
large, near 180°, even if they have similar length. When the 
main peak at around 90° in Cryptomeria was separated into 
two peaks, the distance between them was only 4° (Fig. 5a). 
Although such evaluations will be reported in a following 
report, the large peak was calculated as one floor having 
some angular fluctuations in this report. We therefore cat- 




Fig. 8a-c. Reconstruction of Agathis sp. tracheid cross shape, a, b 
Pentagonal (house) model, and c hexagonal model. Cell side directions 
and standard deviations o are added. Values in parentheses are calcu¬ 
lated from models; b is categorized as a house model although the floor 
contains two sides. House model: hexagonal model = 40.1: 59.9 


egorized them to the house model. The house model was 
reconstructed by the rotation of r : and r 2 of the block model 
around the middle point of each roof side and stopped in 
the position at 59° and 122°. Each side is elongated to the 
crossing points of rotated sides. The hexagonal cell shapes 
are also abundant in the original picture judging from the 
peak area ratio between r : and r 2 and r f , resulting in the 
figure consisting of every two-faced cell wall as a whole. 
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Table 2. Angles of cell walls calculated from the angular distribution function of net map-power spectrum patterns (NM-PSPs) 


Peak angle (0) and o of NM-PSP’s angular distribution function (°) Slope of radial walls 


Agathis sp. 

0 ± 7 

59 ± 12 


91 ± 9 


122 ± 12 

32 

0 

31 

Peak name 

t 

r 2 


r f 


r i 

r 2 

r f 

r i 

Area ratio 

3.1 

1.4 


1 


1.3 




Cryptomeria japonica 

0 ± 5 

76 ± 12 


89 ± 5 


102 ± 12 

13 

0 

13 

Area ratio 

1.9 

1.7 


1 


1.8 




Magnolia obovata 

0 ± 5 

50 ± 7 

70 ± 6 

90 ± 6 

107 ± 6 

131 ± 7 

40 

0 

41 

Peak name 

Peak 1 

Peak 2 

Peak 3 

Peak 4 

Peak 5 

Peak 6 

Peak 2 

r f 

Peak 6 

Area ratio 

0.1 

0.5 

0.3 

1 

0.3 

0.4 





Each power was generated from tangential wall (t) and radial walls (r 2 , r f , and r : ) 


P ac (0) of Magnolia (Fig. 4) was quite different from that 
of Agathis (Fig. 3) and Cryptomeria (Fig. 5a) when it was 
observed. Four peaks are very sharp in the range 0° ^ 0 < 
180°, and small shoulders are also discernable around 75° 
and 105° (arrows in Fig. 4b). Then, it was separated into six 
peaks (peaks 1-6 in Fig. 4c and Table 2) using Sigma Plot, 
suggesting a dominant octagonal cell shape. Cell arrange¬ 
ment and shape are more complicated in hardwoods than in 
softwoods for three main reasons: (1) transverse expansion 
of the vessels, (2) elongation of wood fiber tips, and (3) ray 
distribution. Vessels are very abundant; more than 100/mm 2 
in Magnolia. The cell arrangement and shape are also af¬ 
fected by tip invasion into the cell corner from upper and 
lower fibers. The neighboring radial walls with rays are flat. 
Small invading tips were often observed at the corner posi¬ 
tion between ray and two fibers, and these may have de¬ 
rived the small peaks 3 and 5. In particular, peak 4 had small 
angular fluctuations in the radial walls that were flattened 
by the neighboring rays. When Fig. 4b,c is observed, powers 
dispersed in all directions are equal to the amount of peaks 
1-6, suggesting that many cell sides are not yet identified. 
For the reasons described above, it is difficult to reconstruct 
the correct cell shape of fiber bodies and so only data is 
given in this report. 

Existence ratio of transverse cell shape models 

The peak areas (A peak ) correspond to the accumulation 
length of each cell wall as stated above. Then, we examined 
the existence ratio of transverse cell shape models from 
peak areas in Agathis. We focused on the r wall only, and 
computed the following process in order to simplify calcula¬ 
tions. The area of each peak r 1? r f , r 2 , namely, Ar 1? Ar f , Ar 2 , 
was calculated from the angular distribution function 
(Fig. 3c) after their peak separation. When the existence 
ratio of the house model and hexagonal model were set as 
x%, y %, and wall length in house model as Lr 1? Lr f , Lr 2 , 
respectively, the following relations are: 

x + y = 100 (2) 

(Aii + Ar 2 )/Ar f = (x + ly)^ + Lr 2 )/xLr t (3) 

In Eq. 3, the left side expresses peak area ratio of oblique 
walls (r 1? r 2 ) against the flat wall (r f ) , and the right side 
shows the wall length ratio. From Eqs. 2 and 3 x and y were 


59.9% and 40.1%, respectively (Fig. 8). The fact that the 
house model is much abundant than hexagonal model is 
very attractive. Tracheids with a house model shape existed 
at various places, in addition to both sides of the rays when 
the NM of a transverse section was observed. On the other 
hand, the house model was less abundant than the hexago¬ 
nal model, 45% and 55%, respectively, in Cryptomeria (Fig. 
5a,b). The existence of the rectangular model was ignored 
in the above-mentioned operation process, because the in¬ 
fluence will be only slight because of its rare existence. 

Fluctuation analysis of tracheid shape 

Fluctuations in cell arrangement were analyzed previously 6 . 
Here, the fluctuation in cell wall direction is analyzed. We 
regarded B(x , y), representing the distribution of DM, as 
the periodic function. 6 P(k , 0) was finally translated into the 
periodic distribution function Period (r, 0). Disorders of 
periodicity were expressed as a standard deviation by con¬ 
sidering a normal distribution of the functions. Objects are 
considered as nonperiodic in the direction analysis of the 
cell wall. Each cell side causes two power radii perpendicu¬ 
lar to the cell side from the origin of the PSP. Thus, it was 
sufficient to read P(0). 

We required a position equivalent to the standard devia¬ 
tion o and evaluated the fluctuation in Agathis shown in 
Fig. 3c and Cryptomeria shown in Fig. 5a (Table 2) using a 
normal distribution analysis. 6 Consequently, radial walls r 2 , 
r f , and r t had fluctuations of 12°, 9°, and 12° for sides from 
59°, 91°, and 122°, respectively, while tangential wall t was 
7° smaller than the radial walls in Agathis. On the other 
hand, radial walls are 12°, 5°, and 12° from 76°, 89°, and 
102°, and the tangential wall was 5°, which was equal to r f in 
Cryptomeria. 

Comparison between cell shape models by FTIA 
and OAIA 

Cell arrangement was previously evaluated quantitatively 
by FTIA 6 and cell shape is characterized by FTIA in this 
report. FTIA itself was a new analysis method and there are 
many preconditions. Then, the provability of reconstruction 
of cell shape models should be inspected by OAIA, a con¬ 
ventional analysis method, also using the data of Table 2 in 
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Table 3. Comparison between Fourier transform image analysis (FTIA) and ordinary accumulation image analysis (OAIA) in Agathis sp. 


Parameters 

FTIA 

House model 

Average 

Hexagonal model 

OAIA 

Feret’s H ± o (/mi) 

_ 


_ 

38.9 a ± 6.5 

Feret’s V (/mi) 

— 


— 

49.0 

r-dia compensated from Feret’s H (ptm) 

— 


— 

36.3 b 

r-dia ± o in block model (wm) 

Area (wm 2 ) 

35.l a ± 8.0 

1.46 x 10 3c 

35.l a ± 8.0 

1.51 x 10 3d 

t-dia ± a in block model (/mi) 

41.5 a ± 9.5 


41.5 a ± 9.5 

— 

t-dia minimum («m) 

(35.7, 36.5) 

32.9 

(30.8) 

25 T 

t-dia maximum (jum) 

(46.9) 

50.1 

(52.3) 

— 

h ± a O’) 

31 ± 12 


31 ± 12 

35 ± 9 f 

r 2 ± o (°) 

32 ± 12 


32 ± 12 

38 ± 9 f 

Ratio between two models (%) 8 

40.1 


59.9 

— 


Figures in parentheses are calculated from models 
o, Standard deviation 

a Median value of Gaussian distribution (the others are most frequent values) 

b Compensated value from overestimated Feret’s H (2.6/rni) generated by t-wall inclination (average 84.3°) 
c Calculated with consideration of ratio between two models 
d Compensated value by the addition of half value of IM area 
e Arithmetic mean value 
f Calculated from length sum 

8 Calculated from area ratio of inclining radial walls against flat radial wall at angular distribution function and wall length ratio of models 


our previous report. 6 Agathis was analyzed as an example. 
A very wide region containing about 6600 tracheids was 
used 6 as sample images in FTIA analysis. About 1000 
tracheids were sampled on one screen and then two screens 
were overlapped in this study. On the other hand, 
OAIA measurement became possible only after strict se¬ 
paration of each element. These separations are very 
troublesome. Prudent sampling under higher magnification 
is necessary for the separation. As already described, espe¬ 
cially in OAIA, sampling regions were limited and the 
number of tracheids was about 350. In the pure earlywood 
region, the OAIA sample consists of slightly larger trache¬ 
ids along the radial direction than the tracheids used for 
the FTIA analysis. When this fact is taken into consider¬ 
ation, it can be regarded that the results of both analyses in 
Table 3 are consistent. Some items will be examined in 
more detail. 

In comparison of radial diameters, because tracheids are 
produced constantly from the cambium forming the radial 
files, the radial diameter measured by OAIA is expected to 
be identical to the radial periodicity estimated by FTIA; 
that is, the radial interval of the block model. Feret’s diam¬ 
eter H seems to be the most easily measured (Feret’s H in 
Table 3). Tangential walls are sometimes inclined, and their 
inclinations result in an overestimation of radial diameters. 
Although direction measurements of short lines are very 
difficult in OAIA, especially under low magnification, the 
results for t, r l5 r f , and r 2 are shown in Table 3. Feret’s H data 
were compensated (r-dia compensated from Feret’s H in 
Table 3). As a result, the OAIA measurement of 36.3/un 
was 3.4% larger than the FTIA measurement of 35.1/un. 
The difference must reflect the sample characteristics, be¬ 
cause large tracheids were selectively sampled from pure 
earlywood in OAIA. 

In comparison of cell areas, we examined the validity of 
models from the area information. When transverse sec¬ 


tions of the house and hexagonal models were computed 
considering the existing rates of 59.9% and 40.1%, respec¬ 
tively. Total area average was calculated to be 1.46 X 
10 3 //m 2 . On the other hand, in the OAIA measurement, the 
expanded lumen map (e-LM) 4 separated by intercellular 
layer map (IM) was used. Data must be compensated for by 
the addition of a half value of the IM area surrounding each 
tracheid. Then, it became 1.51 X 10 3 //m 2 and was 3.4% 
larger than FTIA models of 1.46 X 10 3 //m 2 (Table 3). Rays 
are counted in the t-diameter by the periodicity of FTIA, 
and are thought to have weak influence because they oc¬ 
cupy a large part of the xylem. The 3.4% inconsistency 
seems to be due to the reason described above. 

In comparison of cell side lengths, the result was consid¬ 
erably different between OAIA and FTIA (t-dia minimum 
in Table 3). The length of the t-fragment measured by 
OAIA was estimated to be 25.7/un. On the other hand, 
tangential wall length was calculated to be 32.9/un when the 
ratio of both models was considered by FTIA. The result 
of OAIA was smaller by 25%, when both data were 
compared. Very short t-fragments of tracheid tips must be 
the main reason for this underestimation. Moreover, 
fluctuations of x 1 and r 2 in FTIA were 12°. If the angle of r, 
and r 2 is corrected to a slightly steep slope, a tangential 
wall suddenly becomes short. OAIA angles of and r 2 
were larger by 4°-6° according to FTIA. Feret’s V, however, 
was consistent with the t-dia maximum. Determining the 
lengths of the cell sides calculated from angles still has some 
problems. Comparison of the FTIA with OAIA results 
(Table 3) indicates that the data are similar except for only 
a few cases, in spite of the complete independence of both 
analytical processes. Thus, it can be concluded that FTIA is 
appropriate for the quantitative evaluation of softwood tra¬ 
cheids. A large number of structure elements can be ana¬ 
lyzed in a short time by FTIA, so it is suitable for the 
analysis of the most frequent value and also the fluctuation. 
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Fig. 9. General series of FTIA for transverse shape analysis. The pro¬ 
cedure shows maps, patterns and models (in rectangles ), and terms and 
functions from previous reports 4 ' 6 and the present study 


very difficult in hardwoods, however, because rays, vessels, 
and tips disturb cell arrangements and E spn is inevitably 
generated. Some problems remain when the method is used 
for all species, but these problems are expected to be re¬ 
solved. The FTIA series from this report and our previous 
work 6 are shown in Fig. 9 as a conclusion. 
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Then tracheid cell-shape models of many species from all 
over the world will be represented by FTIA in a following 
report. This method can also be applied to other softwoods 
because softwoods consist of very uniform cells. Analysis is 


































